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Abstract 

The recent status of hard thermal photon production in relativistic heavy 
ion collisions is reviewed and the current rates are presented with emphasis 
on corrected bremsstrahlung processes in the quark-gluon plasma (QGP) and 
quark-hadron duality. Employing Bjorken hydrodynamics with an EOS sup- 
porting the phase transition from QGP to hot hadron gas (HHG), thermal 
photon spectra are computed. For SPS 158 GeV Pb-|-Pb collisions, compari- 
son with other theoretical results and the WA98 direct photon data indicates 
significant contributions due to prompt photons. Extrapolating the presented 
approach to RHIC and LHC experiments, predictions of the thermal photon 
spectrum show a QGP outshining the HHG in the high-pT-region. 
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1 Introduction 



Hard real photons are, as dileptons, promising for providing a signature of a quark- 
gluon plasma (QGP) possibly produced in relativistic heavy-ion collisions Direct 
photon spectra have been measured by the WA80 (upper limit) and the WA98 
collaboration at the SPS [21 E] and further data is expected from the RHIC and LHC 
heavy- ion experiments [HIS]. For theoretical investigations on the SPS direct photon 
data and on RHIC and LHC predictions, prompt and thermal photon production 
rates must be convoluted with the space-time evolution of the fireball, while photons 
from the decay of hadrons after freeze-out are already subtracted in the experimental 
analysis. 

In this work, the present status of hard thermal photon production in a QGP 
and a hot hadron gas (HHG) is reviewed critically. Contributions from the QGP 
are presented with emphasis on bremsstrahlung processes that are illustrated for 
the first time in their corrected form. Next, the photon producing processes in the 
HHG are discussed and a conservative estimate of the thermal rate is given. The 
review of the rates is completed by addressing the question of quark-hadron duality. 
With the obtained insights, the recent estimates of the rates are employed together 
with the well-understood Bjorken hydrodynamics [6l [7] to extract the essential fea- 
tures of the thermal and prompt photon yield in the phase transition scenario [8]. 
Comparison with other theoretical work demonstrates the competence of the simple 
hydrodynamical model at SPS energies and substantiates subsequent inspections of 
the WA98 data on direct photon production in SPS 158 GeV Pb+Pb collisions. For 
RHIC and LHC experiments, the rate estimations with the same hydrodynamical 
model are used to predict the pr-range in which the QGP outshines the HHG. 



2 Thermal Photon Production in the QGP 

The production rate for hard (E ^ T) thermal photon^ from an equilibrated QGP 
has been calculated in perturbative thermal QCD applying the hard thermal loop 
(HTL) resummation [10] to account for medium effects. The Compton scattering 
and qq- annihilation contribution 



l—loop 



0.0281 aaMi^^) e"^/^, 



^In the experimentally interesting case E ^ T analytic expressions can be derived since e.g. 
Boltzmann distributions can be used for the outgoing partons [Sj . 
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is derived from the 1-loop HTL photon-polarisation tensor [9l [TT|, [12], while the 
contributions from bremsstrahlung, 

= 0.0219 aa.T^e"-^/^, (2) 

bremss 

and qq- annihilation with an additional scattering in the medium, 

= 0.0105 aasETe-^/^, (3) 

qq—aws 

are obtained from the 2- loop HTL photon-polarisation tensor jl3j, where all three 
rates are listed for a two-flavored {Nf = 2) QGP. Surprisingly, the 2-loop rates, 
([2]) and ([3]), show up at order aas and enhance the spectrum from the QGP phase 
by about a factor of 3 in the experimentally relevant pT-range. Since earlier in- 
vestigations on the importance of bremsstrahlung processes in the QGP \TE\ [TB] 
employed the rates ([2]) and multiplied erroneously by a factor of 4, we present this 
behavior for the flrst time in its corrected formj^ The enhancement is anticipated by 
comparing the rates for flxed temperatures as shown in Fig. 1 and is documented in 
Fig. 2, where the thermal photon spectrum from the QGP, calculated in the model 
discussed below, is presentedjl Traced back to strong collinear singularities [13], this 
2-loop enhancement substantially increased the interest in higher loop contributions. 
In fact, for real photons the 3-loop contribution turns out to be of the same order 
in as as the 2-loop contribution [18]. Since this is very likely the case for higher 
loop contributions as well, one can conclude that thermal photon production in the 
QGP is a non-perturbative mechanism that cannot be accessed in perturbative HTL 
resummed thermal fleld theory. Interestingly, ascribing the quarks a flnite mean 
free path in the QGP, which simulates the Landau-Pomeranchuk-Migdal (LPM) ef- 
fect, does not eliminate important non-perturbative aspects but helps to disarm the 
collinear singularities [19]. This points at possible destructive interferences via the 
LPM effect that do not allow one to interpret the sum of the rates ([I]), ([2]), and ([3]) 
as a lower limit. Adding further uncertainties such as the <^ 1 assumption in the 
HTL calculations contrasted with realistic values of the strong coupling, = 2 — 3, 
one must consider the sum of the rates ([1]), ([2]), and ([3]) only an educated guess. 
However, it is employed in this work as the best result available. Extensions to a 
QGP not in chemical equilibrium, which seems realistic at RHIC and LHC, can also 
be found [201 EB [221 ESI [21] • 

^Due to a miscalculation of the two iV/-dependent constants Jt and Jl (exactly a factor 4 too 
large) in [T31[I1], the derivation of the rates led to Eqs. ^ and ^ both multiplied by a factor 
of 4. In this way, an erroneous enhancement of the spectrum from the QGP phase by about one 
order of magnitude was found in the experimentally relevant pr-range }141 [TSl ITB] . 

•^In the QGP rates, as{T) = (33_2jVj,^Tn(8T/T ) [^ is applied, where the number of flavors 
present in the QGP is set to TV/ = 2 and the transition temperature is set to Tc = 170 MeV. 
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Hard Thermal Photon Rates in the QGP 




E [GeV] 

Figure 1: Hard Thermal Photon Rates in the QGP. For two fixed temperatures, 
T = 150 MeV (lower lines) and T = 250 MeV (upper lines), the QGP hard thermal 
photon rates ([I]), ([2]), and ([3]) are illustrated in the dotted, dashed and dot-dashed 
lines respectively, where the Q;s(r)-parameterization of Karsch [17] is applied. The 
sum indicating the total QGP contribution up to 2-loop order is displayed in the 
solid lines. Inclusion of the 2-loop processes enhances the total QGP rate by about 
a factor of 3. 
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gh = 3 gh = 8 




[GeV] _ [GeV] 



Figure 2: Hard Thermal Photon Spectrum from the QGP. The QGP contributions 
to the hard thermal photon yield from the rates ([T]), ([2]), ([3]), and the sum of these 
are illustrated in the dotted, dashed, dot-dashed and solid lines respectively. The 
integration of the rates over the space-time evolution of the fireball was performed 
in Bjorken hydrodynamics described in the text. Results are shown for two different 
initial temperatures. To = 200 MeV (lower lines) and Tq = 300 MeV (upper lines), 
and two different hadron gas EOS's characterized hy gh = S (left plot) and gh = 8 
(right plot) with the remaining parameters set to the displayed typical values. The 
2-loop processes enhance the total QGP spectrum by about a factor of 3. 
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3 Thermal Photon Production in the HHG 



The thermal photon production in an equihbrated hot hadron gas (HHG) is com- 
puted in effective theories with hadronic couphngs inferred from experiment. A first 
investigation based on a model Lagrangian describing the interaction of vr, p, and 
f] mesons with photons found dominant contributions from the reactions vrvr frf 
and vrp — >■ 7r7 and from the decays a; — >■ 717 and p nir'y [9], [25]. By considering 
additionally the np ^ ai ^ n'j reaction, a strong enhancement of the rate was 
observed [26j. This triggered a more complete and consistent computation using 
an effective chiral Lagrangian with vr, p, and ai mesons. The result was an even 
higher rate due to the presence of the ai meson also in the vrvr — > pj and p —>■ nn'j 
reactions [271 128] ■ Keeping the strong dependence on the model Lagrangian and 
the uncertainties associated with medium effects in mind, the result parameterized 
in [28] supplemented by the uj ^ TC'-f decay |9j can be considered as a conservative 
expression for the emissivity of the HHG. It is employed in this work using the 
formulas listed in the Appendix. For hard photons, E > 1 GeV, we found a rough 
estimate of this sum by multiplying the parameterization given in Eq. (18) of [26] 
by a factor of two 



E 



dN 
d'^x d^p 



4.8 T^-^^ g-l/(1.35T£;)0-77 ^-E/T 



(4) 



had 



where photon energy E and temperature T are to be given in GeV to obtain the 
rate in units of fm~^GeV~^. Results from other Lagrangians can be found in [29|[30] 
and medium effects are investigated in [28l [301 [SD 132], where a significant increase 
of the static rate is observed by dropping the in-medium meson masses that is, 
however, in part compensated in the spectrum through a consequent modification 
of the fireball evolution. Further, the influence of finite chemical potential [31] , finite 
baryon density [33] and additional reactions involving oi, 61, Ki, and other strange 
mesons [Ml [35] have been considered but are neglected in this work. 



4 Quark-Hadron Duality 



In order to test the hypothesis of quark-hadron duality in photon production [361 
[57] . we compare the discussed QGP and HHG hard thermal photon rates for two 
fixed temperatures, T = 150 and 200 MeV, as shown in Fig. 3, where the solid 
and dashed lines indicate the QGP and the HHG contributions respectively. For 
temperatures between 150 MeV and 200 MeV, an energy window can be found, 
where the steeper QGP and the flatter HHG emissivities cross. Considering the 
uncertainties in the rates, the validity of the classic statement 'the hadron gas shines 
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Figure 3: Quark-Hadron Duality in the Hard Thermal Photon Rates. The sum of 
the QGP contributions ([1]), ([2]), and ([3]), illustrated in the solid lines, and that of 
the HHG contributions (Song's results [271 [2H] supplemented by the u de- 
cay [9]), illustrated in the dashed lines, are compared for two fixed temperatures, 
T = 150 MeV (lower lines) and T = 200 MeV (upper lines), where the a<j(T)- 
parameterization of Karsch [TJ| is applied in the QGP rates. Considering the un- 
certainties in the rates, quark-hadron duality in the rates at a certain temperature 
can neither be confirmed nor be ruled out. 
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as brightly as the quark-gluon plasma' [9] can neither be confirmed nor be ruled out, 
at least in a certain temperature regime, by taking into account the 2-loop QGP 
processes and the Oi meson in the HHG. However, as indicated by Fig. 3, there is 
no reason to assume that the photon production rates from the two phases coincide 
at all temperatures. But even if the rates of the two phases agree in the relevant 
temperature regime, the spectrum might be changed by the presence of the QGP 
due to a different space-time evolution. 



5 Thermal Photon Spectrum 

The observable quantity is the spectrum of direct photons which is calculated the- 
oretically by convoluting the prompt and thermal photon rates with the space-time 
evolution of the fireball produced in heavy-ion collisions. Concentrating on the 
effects of the thermal photon rates, we describe the thermalized collision phase in 
Bjorken-l+l-hydrodynamics [SIE] taking an EOS with a first-order phase transition 
from a QGP, modelled by a massless two-flavored parton gas, to a HHG, modelled by 
a massless pion gas. This well-understood, simple approach allows us to demonstrate 
very clearly the influences of the rates on the spectrum, and interestingly, compar- 
ison with other work shows it to be as competent as 2+1 hydrodynamical models 
at least at SPS energies pS]. This means that at SPS the transverse expansion has 
only minor effects on the photon production. 

5.1 Comparison with other Work 

The WA98 direct photon data analysis of Gallmeister et al. obtained in a model 
describing a spherically symmetric expansion [37] can be reproduced in the simple 
model with the discussed ratesQ By assuming a thermalization time of tq = 1 fm and 
equal initial, transition, and freeze-out temperatures of Tq = Tc = Tf = 170 MeV 
as in [37], i.e. only a mixed phase, this single temperature value suffices to describe 
the WA98 data [21 E] when the prompt photon estimation of [37] is added. 

It is also possible to reproduce the WA98 direct photon data analysis of Sri- 
vastava et al. ^39j, which does not necessitate prompt photons but instead initial 
conditions that are rather extreme for SPS, i.e. a very small thermalization time 
of To = 0.2 fm and a very high initial temperature of Tq = 335 MeV. Since their 
analysis was performed using the 2-loop QGP rates by mistake multiplied by a fac- 
tor of 4, even smaller thermalization times together with higher initial temperatures 
are demanded when applying the corrected 2-loop QGP rates. This is due to the 

^The uj TTj decay is not taken into account since electromagnetic meson decays were sub- 
tracted in the experimental analysis of the WA98 collaboration j38j . 
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relative importance of the QGP rates in the high p^-range. Concentrating on a 
vahdity-check of the simple model, a comparison using the same erroneous rates is 
performed. While the QGP thermal photon spectrum obtained in our simple model 
matches directly the one of Srivastava et al. we must increase the effective de- 
grees of freedom in the ideal massless pion gas from the actual value of Qh = 3 to an 
effective one of gh = 8 in order to achieve the fit in the HHG thermal photon spec- 
trum. This points to the rich HHG EOS employed in [39j having a much stronger 
effect than transverse expansion. 



5.2 Thermal Photon Spectrum from the QGP 

These insights into the reliability of the simple model substantiate the QGP thermal 
photon spectrum calculated in this model and shown in Fig. 2 to illustrate the effect 
of 2-loop processes in the QGP. Displaying the spectrum not only for gh = 3 (left 
plot) but also for gt = 8 (right plot) demonstrates additionally the effect of a richer 
HHG EOS, which lowers the QGP thermal photon yield by reducing the lifetime of 
the mixed phase. 



5.3 Comparison with WA98 Direct Photon Data 

The total thermal photon spectrum emerging from the discussed rates in the simple 
model can be exploited by using the WA98 measurements of the direct photon 
yield [3] to specify upper limits on the initial temperature, T^""^ , reached in the 
SPS 158 GeV Pb-|-Pb collisions. This is presented in Fig. 4. For typical parameters 
To = 1 fm, Tc = 170 MeV, Tj = 150 MeV, nucleon number A = 208 (corresponding 
to Pb + Pb collisions), projectile rapidity ynud = 2.8 (corresponding to a center- 
of-mass energy of ^/s = 17 GeV), and with an ideal massless pion gas of (^/i = 3 
effective degrees of freedom, even with an initial temperature of Tq = Tc = 170 MeV 
(dotted line) the computed spectrum exceeds the experimental limits at low ptI 
Simulating the richer HHG EOS by setting gh = 8, the lifetime of the mixed and 
pure HHG phase drops, resulting in a decrease of the spectrum that brings back 
the possibility of having a phase transition scenario with T^"^ = 220 MeV (short 
dashed line). By reducing the critical temperature to = 160 MeV, the mean 
temperature is reduced enabling also for gh = 3 the phase transition scenario with 
rp^ax ^ MeV (solid line). For gh = 8 and = 160 MeV, T^^"^ = 235 MeV 
(long dashed line) is found, which might be interpreted as an upper bound on the 
initial temperature reached in the SPS 158 GeV Pb+Pb colhsionsEl However, the 

^Due to the mentioned possible destructive interferences via the LPM effect in the QGP, the 
presented J"™"^ values cannot represent rigid upper bounds. They rather indicate the magnitude 
and point to the relative differences caused by different model scenarios. 
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WA98 Direct Photon Data & Thermal Spectrum 
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Figure 4: WA98 Direct Photon Data and Upper Limits on the Initial Temperature. 
The theoretical spectra supporting the phase transition from QGP to HHG (solid, 
short-dashed, and long-dashed lines) are compared with the experimental upper 
limits (vertical lines) and data (points with error bars) . The discrepancy in the high- 
PT-region demonstrates that a significant prompt photon contribution is necessary 
for a theoretical explanation of the experimental results within the phase transition 



scenario. 
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1.0 


190 
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5.3 
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0.5 


310 


LHC 


ALICE 


208 


5500 


8.6 


5625 


0.5 


450 



Table 1: Parameters used in the Predictions of the Thermal Photon Spectrum 
for RHIC and LHC Experiments. The listed values for multiplicity dN/dy, 
thermalization time Tq, and initial temperature Tq are the estimates given 
in [H]. 

spectra supporting the phase transition from QGP to HHG cannot describe the 
direct photon data for pt > 2 GeV as shown in Fig. 4. This discrepancy in the high- 
PT-region demonstrates that a significant prompt photon contribution is necessary 
for a theoretical explanation of the experimental results within the phase transition 
scenario. 

5.4 Predictions for RHIC and LHC 

Finally, we employ the simple model to extract predictions for RHIC and LHC. 
Although earlier comparisons with |T3] revealed a more important transverse ex- 
pansion [TB], it is interesting to locate in the above approach the p^-range in which 
the QGP outshines the HHG. Assuming a three-flavored QGP, the prefactors on 
the rhs in Eqs. (P), (ED, and become 0.0338, 0.0281, and 0.0135 respectively. 
Further, the SPS, RHIC, and LHC parameters of [H] are adopted as summarized 
in Tab. 1. Figure 5 shows our results for the QGP, HHG, and sum of both contri- 
butions in the dashed, dotted, and solid lines respectively. The spectra for gn = ^ 
are shown on the left and the ones for gh = 8 on the right. As a reference point, 
the SPS result obtained with the parameters of [2] is also presented. For (7/1 = 3 
{gh = 8), the QGP outshines the HHG for pr > 3.2 GeV {pr > 2.2 GeV) at RHIC 
and for pr > 2.5 GeV {pr > 1.7 GeV) at LHC, while at SPS photons from HHG 
dominate the spectrum for gh = ^ and gh = 8 aX all pt's. Since the employed simple 
model does not describe deviations from chemical equilibrium, the reader is referred 
to [21] , where this issue is addressed systematically and RHIC and LHC predictions 
are provided for the QGP photon emissivity. 
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gh = 3 g, = 8 




Pt [GeV] - [GeV] 



Figure 5: QGP vs. HHG Photon Emissivities in SPS, RHIC, and LHC Pb+Pb 

Collisions. For gn = S {qh = 8), the QGP outshines the HHG for pr > 3.2 GeV 
{Pt > 2.2 GeV) at RHIC and for pr > 2.5 GeV {pr > 1.7 GeV) at LHC, while 
at SPS photons from HHG dominate the spectrum for (7/1 = 3 and gh = 8 at all 
PtS. The initial conditions vary from (tq = 1 fm, Tq = 190 MeV) at SPS over 
(to = 0.5 fm, To = 310 MeV) at RHIC to (tq = 0.5 fm, Tq = 450 MeV) at LHC, 
while the same values of — 160 MeV and Tf — 100 MeV are chosen. 
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6 Conclusion 



We have discussed the hard thermal photon production in relativistic heavy ion 
colhsions using the most recent estimates for the photon production rates from 
the QGP and the HHG. For the first time, the enhancement of the QGP thermal 
photon rate and yield due to bremsstrahlung processes was presented in its corrected 
form. The consideration of these processes hinted at the recently confirmed non- 
perturbative nature of thermal photon production in the QGP [18]. The derivation 
of a reliable rate requires now new developments in thermal field theory [19]. For 
thermal photon production in the HHG, the strong dependence on assumptions 
regarding the thermalized hadron species, their interactions, and the role of medium 
effects was emphasized. We identified the best estimates available for the QGP and 
HHG rates and presented a convenient parameterization for the total HHG rate. 
Using these estimates, we could find no indication of a quark-hadron duality in the 
photon production rate. 

Integrating the rate estimations over the space-time evolution of the fireball 
modelled in Bjorken hydrodynamics, we obtained the thermal photon spectrum. 
Comparison with other models, which describe transverse expansion of the fireball, 
demonstrated the competence of the simple model for the phase transition scenario 
at SPS energies and substantiated the subsequent investigation of the WA98 direct 
photon data. We found that this experimental data allows no conclusion to be drawn 
about the existence of a QGP phase at SPS but can be explained by a conservative 
thermal source plus prompt photons. Finally, the simple model was employed also 
for RHIC and LHC energies and predicted a QGP outshining the HHG in the high- 
PT-range. If this picture is not spoiled by transverse expansion and chemical non- 
equilibrium or covered by dominant prompt photon contributions, the RHIC and 
LHC experiments might see thermal photons from the QGP in a certain photon 
momentum range that could provide the desired signature of the QGP. 

As RHIC is already taking data and LHC under construction, a reliable expres- 
sion for the thermal photon production rate in the QGP is of utmost importance. 
Optimistically the final result could reveal photons to be a smoking gun for the 
production of the QGP. In order to reduce the uncertainties concerning the fireball 
evolution, we recommend a combined investigation of real photon, dilepton, and 
hadron spectra, similar to the one of Sollfrank et al. [ID], in which aspects such 
as different EOS's, finite baryon density, chemical non-equilibrium, and transverse 
expansion should be addressed in a systematic way. Once such an investigation is 
completed in accordance with the SPS data, it should be extended to RHIC and 
LHC energies to obtain serious predictions. 
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Appendix 



The employed rate for thermal photon production in the HHG is composed of the 
exact expression for the decay — 7r7 ^ and the parameterizations for the processes 
TTTT fTf, up — i> 7r7, and p n-K'y, in which the ai meson is taken into account 
properly [271 12H] • For completeness and convenience, the explicit formulas are listed. 
By inserting E, T, and the pion mass m^r in GeV, the following parameterizations 
for process = vrvr — p7, Tip — >■ 717, and p — > inc'-f reproduce the corresponding rates 
in units of fm~^GeV~^ 



d^x d^p 



T^e ^'"^ Fprocess{T/m^,E/m^), (5) 



process 



where 



y) = exp[-12.055 + 4.387x + (0.3755 + 0.00826x)?/ 

+ (-0.00777 + 0.000279a;)?/2 + (5.7869 - 1.0258x)/?/ 

+ (-1.979 + 0.58x)/?/2], (6) 

F^p^^^ix, y) = exp[-2.447 + 0.796x + (0.0338 + 0.0528x)y 

+ (-21.447 + 8.2179x)/y + (1.52436 - 0.38562a;) /y^], (7) 

Fp^^^^ix, y) = exp[-6.295 + 1.6459x + (-0.4015 + 0.089x)?/ 

+ (-0.954 + 2.05777x)/y]. (8) 



No parameterization is used for the decay of the u meson since the exact expression 
requires only a one-dimensional integration [9] 



^.93x10"^ GeV 



E 



/ dE^E^fB{E^)[l + fB{E^-E)], (9) 



d^x d^p 

where E^in = 1-03 {E^ + 0.14 GeV^)/E and /b(^) = l/[exp(E/T) - 1]. 
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Figure 2: Hard Thermal Photon Spectrum from the QGP. The QGP contributions 
to the hard thermal photon yield from the rates (1), (2), (3), and the sum of these 
are illustrated in the dotted, dashed, dot-dashed and solid lines respectively. The 
integration of the rates over the space-time evolution of the fireball was performed 
in Bjorken hydrodynamics described in the text. Results are shown for two different 
initial temperatures, Tq = 200 MeV (lower lines) and Tq = 300 MeV (upper lines), 
and two different hadron gas EOS's characterized by gh = 3 (left plot) and gh = 8 
(right plot) with the remaining parameters set to the displayed typical values. The 
2-loop processes enhance the total QGP spectrum by about a factor of 3. 
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Figure 4: WA98 Direct Photon Data and Upper Limits on the Initial Temperature. 
The theoretical spectra supporting the phase transition from QGP to HHG (solid, 
short-dashed, and long-dashed lines) are compared with the experimental upper 
limits (vertical lines) and data (points with error bars) . The discrepancy in the high- 
PT-region demonstrates that a significant prompt photon contribution is necessary 
for a theoretical explanation of the experimental results within the phase transition 



scenario. 
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Figure 5: QGP vs. HHG Photon Emissivities in SPS, RHIC, and LHC Pb+Pb 

Collisions. For gn = S {qh = 8), the QGP outshines the HHG for pr > 3.2 GeV 
{Pt > 2.2 GeV) at RHIC and for pr > 2.5 GeV {pr > 1.7 GeV) at LHC, while 
at SPS photons from HHG dominate the spectrum for (7/1 = 3 and gh = 8 at all 
PtS. The initial conditions vary from (tq = 1 fm, Tq = 190 MeV) at SPS over 
(to = 0.5 fm, To = 310 MeV) at RHIC to (tq = 0.5 fm, Tq = 450 MeV) at LHC, 
while the same values of — 160 MeV and Tf — 100 MeV are chosen. 
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